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1 Introduction 



The quantum field theory in curved space-time has been a matter of great interest in 
recent years because of its applications to cosmology and astrophysics. The evidence of 
existence of strong gravitational fields in our Universe led to the study of the quantum 
effects of material fields in external classical gravitational field. After the appearance of 
Parker's paper on scalar fields [[|] and spin-^ fields ||, several authors have studied this 
subject. Although the Universe seems homogenous and isotropic at present, there are 
no observational data guarantying the isotropy in the era prior to the recombination. In 
fact, there are theoretical arguments that sustain the existence of an anisotropic phase 
that approaches an isotropic one ||. Interest in studying Klein-Gordon and Dirac 
equations in anisotropic models has increased since Hu and Parker have shown that 
the creation of scalar particles in anisotropic backgrounds can dissipate the anisotropy 
as the Universe expands. 

A Bianchi type-I (B-I) Universe, being the straightforward generalization of the flat 
Robertson- Walker (RW) Universe, is one of the simplest models of an anisotropic Uni- 
verse that describes a homogenous and spatially flat Universe. Unlike the RW Universe 
which has the same scale factor for each of the three spatial directions, a B-I Universe 
has a different scale factor in each direction, thereby introducing an anisotropy to the 
system. It moreover has the agreeable property that near the singularity it behaves 
like a Kasner Universe even in the presence of matter and consequently falls within the 
general analysis of the singularity given by Belinskii et al ||. And in a Universe filled 
with matter for p = je, 7 < 1, it has been shown that any initial anisotropy in a B-I 
Universe quickly dies away and a B-I Universe eventually evolve into a RW Universe 
||. Since the present-day Universe is surprisingly isotropic, this feature of the B-I 
Universe makes it a prime candidate for studying the possible effects of an anisotropy 
in the early Universe on present-day observations. In light of the importance of men- 
tioned above, several authors have studied linear spinor field equations J/J, || and the 
behavior of gravitational waves (GWs) 0, [10], [1 L| in B-I Universe. Nonlinear spinor 



field (NLSF) in external cosmological gravitation field was first studied by G. N. Shikin 
in 1991 |r2"] . This study was extended by us for more general case where we consider 
nonlinear term as an arbitrary function of all possible invariants generated from spinor 
bilinear forms. In that paper we also studied the possibility of elimination of initial 
singularity specially for Kasner Universe In a recent paper |[14[] we studied the 



behavior of self-consistent NLSF in B-I Universe that was followed by the papers Jl5] 



16H where we studied the self-consistent system of interacting spinor and scalar fields. 



The purpose of the paper is to extend our study for more general NLSF in presence of 
perfect fluid. In the section 2 we derive fundamental equations corresponding to the 
Lagrangian for the self-consistent system of spinor and gravitational fields in presence 
of perfect fluid and seek their general solutions. In section 3 we give a detail analysis 
of the solutions obtained for different kinds of nonlinearity. In section 4 we study the 
role of perfect fluid and in section 5 we sum up the results obtained. 
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2 Fundamental equations and general solutions 



The Lagrangian for the self-consistent system of spinor and gravitation fields in presence 
of perfect fluid is 



R i 
Yk + 2 



r. 



(2.1) 



with R being the scalar curvature and k being the Einstein's gravitational constant. 
The nonlinear term describes the self- interaction of spinor field and can be presented 
as some arbitrary functions of invariants generated from the real bilinear forms of spinor 
field having the form: 

S = P = ^ 7 5 <0, if = (^7 M V), Ail = ($tVV0> = (^cr^if;), 

where a^ v = (i/2)[j ft j v — 7 I/ 7**]. Invariants, corresponding to the bilnear forms, look 

i = s 2 , j = p 2 , i v = v lt v»= $-fi/>) M^/V), 

I A = A^A» = (fyVVO M^VVO, It = T, u = tfa^) v^(^ Q »- 

According to the Pauli-Fierz theorem WJ\ among the five invariants only I and J are 
independent as all other can be expressed by them: I v = —I a = I + J and It = I — J. 
Therefore we choose the nonlinear term Ln = F(I, J), thus claiming that it describes 
the nonlinearity in the most general of its form. L m is the Lagrangian of perfect fluid. 

We choose B-I space-time metric in the form 



ds 2 



dt 2 — jij (t) dx % dx-' . 



(2.2) 



As it admits no rotational matter, the spatial metric jij(t) can be put into diagonal 
form. Now we can rewrite the B-I space-time metric in the form ||18|| : 



ds 2 = dt 2 - a 2 (t) dx 2 - b 2 (t) dy 2 - c 2 (t) dz 2 



(2.3) 



where the velocity of light is taken to be unity. Einstein equations for a(t),b(t) and 
c(t) corresponding to the metric (2.3) and Lagrangian (2.1) read [|TjJ: 
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(2.4) 
(2.5) 
(2.6) 
(2.7) 



where points denote differentiation with respect to t, and T = T£. 
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NLSF equations and components of energy-momentum tensor for the spinor field and 
perfect fluid corresponding to (2.1) are 

H M V M V - rnip + FfiSip + Fj2Pij 5 if; = 0, 

«V^7 M + rmp - FrfSip - Fj2Pi^ 5 = 0, (2.8) 

where Fj := dF/dl and Fj := dF/dJ. 

T£ = \g pv {^l^A + fy„V„V - V^ 7 ^ - V^t^) - 6?L sp + T; (m) , (2.9) 
while L sp on account of spinor field equations takes the form: 

: -[21 F! +2JFj-L N ]. 



. dL N dL N 



Here T^, m ^ is the energy-momentum tensor of perfect fluid. For a Universe filled with 
perfect fluid, in the concomitant system of reference (u° = 1, u l = 0, % = 1,2,3) we 
have 

T; [m) = (p + e)u/ - 5> = (e, -p, -p, -p), (2.10) 

where energy e is related to the pressure p by the equation of state p = 7 e, the general 
solution has been derived by Jacobs |J. 7 varies between the interval < 7 < 1, 
whereas 7 = describes the dust Universe, 7 = | presents radiation Universe, | < 
7 < 1 ascribes hard Universe and 7 = 1 corresponds to the stiff matter. In (2.8) and 
(2.9) V M denotes the covariant derivative of spinor, having the form [IS]: 

V^ = |^-I>, (2.11) 

where T M (x) are spinor affine connection matrices. 7 M (x) matrices are defined for the 
metric (2.3) as follows. Using the equalities pO| , f2~l|] 



9^{x) = e*(x)e b u (x)ri ab , ^(x) = e a ^x)-? a , 

where r] a b = diag(l, — 1, — 1, — 1), 7« are the Dirac matrices of Minkowski space and 
e^(x) are the set of tetradic 4- vectors, we obtain the Dirac matrices 7 M (x) of curved 
space-time 

7 = 7°, 7 1 =7 1 /a(t), l 2 = f/Kt), 7 3 = 7 3 Mt), 
7o = 7o, 7i = 7ia(*)> 72 = l2b(t), 73 = %c{t). 
T^x) matrices are defined by the equality 



r^) = ^(x)(^-r;)7V, 



which gives 



r = o, r\ = \a{t)ff, r 2 = h(t)ff, r 3 = h(t)f^. (2.12) 
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Flat space-time matrices we choose in the form, given in [22 
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Defining 7 s as follows 
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'9 e iivapi ^0123 



-^7°7 1 7 2 7 3 



.^0-1-2-3 = - Bj 



we obtain 
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( 0-10^ 

-1 

-10 

V -1 / 



We study the space-independent solutions to NLSF equation (2.8). In this case the 
first equation of the system (2.8) together with (2.11) and (2.12) is 

%1 \di + ^)^- m ^ + ^ + ^7 5 ^ = 0, r{t) = a(t)b(t)c(t), (2.13) 

where we denote T> := 25* Fj and Q := IP Fj. For the components tp p = V p (t), where 
p = 1, 2, 3, 4, from (2.13) one deduces the following system of equations: 



Vi + 7T V 1 + i 
It 


[m 


-v)Vi-gv 3 


= 0, 


V 2 + ^V 2 + i 
It 


[m 


-v)v 2 - gv 4 


= 0, 


v 3 + ^-V 3 - 1 

It 


[m 


-v)v 3 + gv x 


= 0, 


V a + ^-V 4 - i 
It 


[m 


-v)v 4 + gv 2 


= 0. 



Let us now define the equations for 



p = i(v x v* - v*v 3 + v 2 v: - V*V 4 ), 
r = <y x vz + v*v 3 + v 2 v; + V*V 4 ), 
s = {yzv x + v*v 2 - v;v 3 - v;v 4 ). 



(2.14) 



(2.15) 
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After a little manipulation one finds 



dS 
~dt 



2GR 



dR 
~df 



+ 2(m - V) P + 2QS Q 



dPo 
dt 



2{m-V) R 



0, 
0, 
0, 



(2.16) 



where So = rS, P = tP, Rq = tR. From this system we obtain 

S S Q + RqRo + PoPo = 0, 

that gives 

S 2 + R 2 + P 2 = C 2 /t 2 , C 2 = const. (2.17) 

Let us go back to the system of equations (2.14). It can be written as follows if one 
defines W a = yfr V a : 



Wx + i^Wx-gWz = 0, 
Ws-iQWu + GW! = 0, 



w 2 + i$w 2 - gw 4 = o, 
w 4 - i<s>w A + gw 2 = o, 



(2.18) 



where $ = m — T>. Defining U(a) = W(t), where a — J Qdt, we rewrite the foregoing 
system as: 



ui + i^/g^-Ua = o, u^ + i($/g)u 2 -Ut = o, 
u^-t{^/g)u 3 + u 1 = o, ui-i(Q/g)u A + u 2 = o, 



(2.19) 



where prime (') denotes differentiation with respect to a. 

Let us now solve the Einstein equations. To do it we first write the expressions for 
the components of the energy-momentum tensor explicitly. Using the property of flat 
space-time Dirac matrices and the explicit form of covariant derivative V M one can 
easily find 

T$ = mS- F(I, J) + e, = T| = Tf = 21 Fj + 2 J Fj - F(I, J) - p. (2.20) 
Summation of Einstein equations (2.4), (2.5) and (2.6) leads to the equation 

- = -niTl+Ti + Tl-^T) = ^(mS + 2IF I + 2J Fj - 2F(I, J) +e-p). (2.21) 
r 2 2 

In case if the right hand side of (2.21) be the function of r(t) = a(t)b(t)c(t), this 
equation takes the form 

f + $(r)=0. (2.22) 

As is known this equation possesses exact solutions for arbitrary function $(r). Giv- 
ing the explicit form of Ln = F(I, J), from (2.21) one can find concrete function 
r(t) = abc. Once the value of r is obtained, one can get expressions for components 
V a (t), a = 1,2, 3, 4. Let us express a, b, c through r. For this we notice that subtrac- 
tion of Einstein equations (2.4) - (2.5) leads to the equation 



a 
a 



b dc 
b ac 



be d / d 
be dt\a 



+ 



b c 
bj \a b c 



= 0. 



(2.23) 
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Equation (2.23) possesses the solution 



a ( f dt\ 

— = Diex.pl Xi J — I, D 1 = const., X 1 = const. 



(2.24) 



Subtracting equations (2.4) - (2.6) and (2.5) - (2.6) one finds the equations similar to 
(2.23), having solutions 



(2.25) 



where D 2 , D 3 , X 2 , X 3 are integration constants. There is a functional dependence be- 
tween the constants Di, D 2 , D 3 , X\, X 2) X 3 : 

D 2 = D 1 D 3 , X 2 = X 1 + X 3 . 

Using the equations (2.24) and (2.25), we rewrite a(t) , b(t) , c(i) in the explicit form: 

-2X,+X 3 



a(t) = (DlD^rzexp 
b(t) = (D^ l D 3 )^Aexp 
c(t) = (D 1 D 2 3 )-hhxp 



J 



dt 



3 J r{t) 
X\ — X 3 f dt 



r dt 
J W) 



3 J r(t). 

Xi + 2X 3 r dt 



r dt 

J W) 



(2.26) 



Thus the previous system of Einstein equations is completely integrated. In this process 
of integration only first three of the complete system of Einstein equations have been 
used. General solutions to these three second order equations have been obtained. The 
solutions contain six arbitrary constants: Di, D 3 , X 1 , X 3 and two others, that were 
obtained while solving equation (2.22). Equation (2.7) is the consequence of first three 
of Einstein equations. To verify the correctness of obtained solutions, it is necessary to 
put a,b,c into (2.7). It should lead either to identity or to some additional constraint 
between the constants. Putting a, b, c from (2.26) into (2.7) one can get the following 
equality: 



1 

3r 



+ x l x 3 + x\ 



- K (T -V), (2.27) 



that guaranties the correctness of the solutions obtained. In fact we can rewrite (2.21) 
and (2.27) as 

I = yCC + T, 1 ), (2.28) 



and 



•2 



T 



(2.29) 



where X := X\ + X X X 3 + X\. Combining (2.28) and (2.29) together one gets the 
solution for r in quadrature: 



dr 



^3kt 2 T$ + X/3 



t. 



(2.30) 
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Let us note that in our further study we exploit the equations (2.21) to obtain r and 
(2.27) to estimate integration constants. 

It should be emphasized that we are dealing with cosmological problem and our main 
goal is to investigate the initial and the asymptotic behavior of the field functions and 
the metric ones. As one sees, all these functions are in some functional dependence 
with r: ip ~ an d o« ~ T 1 / 3 e ± I dt / T _ Therefore in our further investigation we 

mainly look for r, though in some particular cases we write down field and metric 
functions explicitly. 

3 Analysis of the solutions obtained for some spe- 
cial choice of nonlinearity 

Let us now study the system for some special choice of L N . First we analyze the system 
only for the NLSF which will be followed by the study when the Universe is filled with 
perfect fluid. But first of all we study the linear case. The reason to get the solution 
to the self-consistent system of equations for the linear spinor and gravitational fields 
is the necessity of comparing this solution with that for the system of equations for the 
nonlinear spinor and gravitational fields that permits to clarify the role of nonlinear 
spinor terms in the evolution of the cosmological model in question. Using the equation 
(2.21) one gets 

r(t) = (1/2) Mt 2 + Vl t + y (3.1) 

where M = ^nmCo, Cq = C\ + Cf — C 2 — C\ and y±,yo are the constants. In this 
case we get explicit expressions for the components of spinor field functions and metric 
functions: 

V r (t) = (C r /V^)e- imt , r = l,2; V t (t) = (C,/v^) e imt , 1 = 3,4- (3.2) 

a(t) = (DlD 3 )k^Mt 2 + yi t + y^Z^ +x ^ B , 

b(t) = (D- 1 D 3 ) 1 s(^Mt 2 + y 1 t + y ) 1 sZ-^- x ^ B , 

c(t) = (D.Dly^^Mf + y.t + yo^Z-^ 2 ^, (3.3) 

where Z = B = M(t 1 - t 2 ), and t lj2 = -yi/M ± ^{y^M) 2 - 2y /M are the 

roots of the quadratic equation Mt 2 + 2yit + 2y = 0. Substituting r(t) into (2.27), 
one gets 

y\ - 2My = (X 2 + X t X 3 + X 3 2 )/3 = X/3 > 0. (3.4) 

This means that the quadratic polynomial in (3.1) possesses real roots, i.e. r(t) in 
(3.1) turns into zero at t — t 1>2 and the solution obtained is the singular one. Let us 
now study the solutions (3.1) - (3.3) at t — > oo. In this case we have 

3 

r(t) w -KmC t 2 , a(t) « b(t) « c(t) w t 2/s , 
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that leads to the conclusion about the asymptotical isotropization of the expansion 
process for the initially anisotropic B-I space. Thus the solution to the self-consistent 
system of equations for the linear spinor and gravitational fields is the singular one 
at the initial time. In the initial state of evolution of the field system the expansion 
process of space is anisotropic, but at t — > oo there happens isotropization of the 
expansion process. 

Once the solutions to the linear spinor field equations and corresponding to them metric 
functions are obtained, let us now study the nonlinear case. 

I. Let us consider the case when = F(I). It is clear that in this case Q = 0. From 
(2.16) we find 

S = C /t, C = const. (3.5) 

As in the considered case L N = F depends only on S, from (3.5) it follows that F(I) 
and Fi(I) are functions of r = abc. Taking this fact into account, integration of the 
system of equations (2.14) leads to the expressions 

V r (t) = (C r /y/¥) e-*\ r = l,2, Vftt) = (Q/y/r) e* n , Z = 3,4. (3.6) 

where C r and C\ are integration constants. Putting (3.6) into (2.15) one gets 

S={Cl + Cl-Cl-Cl)/r. (3.7) 

Comparison of (3.5) with (3.7) gives C = Cf + Cf - Cf - Cf . 

Let us consider the concrete type of NLSF equation with F(I) = \I^ n / 2 ^ = \S n where 
A is the coupling constant, n > 1. In this case for r one gets: 

r = (3/2)«C [m + A(n - 2)C^~ 1 / V 71 " 1 ]. (3.8) 

The first integral of the foregoing equation takes form: 

f 2 = 3kC [mr - A C^/t 71 ' 2 + g% (3.9) 

where from (2.27) one determines g 2 = X/9kCq. The sign Cq is determined by the 
positivity of the energy- density T® of linear spinor field: 

T ° = mC /r>0. (3.10) 

It is obvious from (3.10) that C > 0. Now one can write the solution to the equation 
(3.9) in quadratures: 



/ 



-(n-2)/2 rf7 



mT n-l + g 2 T n-2 _ XQn- 



= = ^kC t (3.11) 



The constant of integration in (3.11) has been taken zero, as it only gives the shift of 
the initial time. Let us study the properties of solution to equation (3.8) for n > 2. 
From (3.11) one gets 

r(t) |^oo~ (3/A)KmC t 2 , (3.12) 

which coincides with the asymptotic solution to the equation (3.3). It leads to the 
conclusion about isotropization of the expansion process of the B-I space. It should 
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be remarked that the isotropization takes place if and only if the spinor field equation 
contains the massive term [cf. the parameter m in (3.12)]. If m=0 the isotropization 
does not take place. In this case from (3.11) we get 



r(t) It^oo* pKC gH. (3.13) 

Substituting (3.13) into (2.26) one comes to the conclusion that the functions a(t), b(t) 
and c{t) are different. Let us consider the properties of solutions to equation (3.8) 
when t — > 0. For A < from (3.11) we get 

r(t) = [(3/4)n\|A|g 1/n t 2/n -> 0, (3.14) 

i.e. solutions are singular. For A > 0, from (3.11) it follows that r = cannot 
be reached for any value of t as in this case the denominator of the integrand in 
(3.11) becomes imaginary. It means that for A > there exist regular solutions to 



the previous system of equations ||14|| . The absence of the initial singularity in the 
considered cosmological solution appears to be consistent with the violation for A > 0, 
of the dominant energy condition in the Hawking- Penrose theorem |18| . 



Let us consider the Heisenberg-Ivanenko equation when in (3.8) n=2 [23|. In this case 
the equation for r(i) does not contain the nonlinear term and its solution coincides 
with that of the linear equation (3.3). With such n chosen the metric functions a, b, c 
are given by the equality (3.2), and the spinor field functions are written as follows: 

V r = (C r /VT) e -^ Z AlXCo/B , Vi = (Q/^) e ^z- AlXCo/B (3.15) 

As in the linear case, the obtained solution is singular at initial time and asymptotically 
isotropic as t — ► oo. 

We now study the properties of solutions to equation (3.8) for 1 < n < 2. In this case 
it is convenient to present the solution (3.11) in the form: 



/ 



dr 



mr - Ar 2 -^" 1 + g 2 



3nC t (3.16) 



As t — > oo, from (3.16) we get the equality (3.12), leading to the isotropization of the 
expansion process. If m = and A > 0, r(t) lies on the interval 

< r(t) < {g 2 /\C^ l ) 1/{2 - n) . 
If m=0 and A < 0, the relation (3.16) at t — > oo leads to the equality: 

r(t) w [(3/4)n 2 K\\\C%} 1/n t 2/n . (3.17) 

Substituting (3.17) into (2.26) and taking into account that at t — > oo 

f dt ^ n(3K\\\n 2 CZ)V" t _ 2/n+1 ^ 
i r~ (n-2)2 2 / n 

due to —2/n + 1 < 0, we obtain 

a(t) ~ b(t) ~ c(t) ~ [r(t)} 1/3 ~ t 2/3n -> oo. (3.18) 
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It means that the solution obtained tends to the isotropic one. In this case the 
isotropization is provided not by the massive parameter, but by the degree n in the 
term L N = XS n . (3.16) implies 



r(t) U_^o~ ^KC o9 n^O, (3.19) 

which means the solution obtained is initially singular. Thus, for 1 < n < 2 there exist 
only singular solutions at initial time. At t — > oo the isotropization of the expansion 
process of B-I space takes place both for m/fl and for m = 0. 

Let us finally study the properties of the solution to the equation (3.8) for < n < 1. 
In this case we use the solution in the form (3.16). As now 2 — n > 1, then with 
the increasing of r(i) in the denominator of the integrand in (3.16) the second term 
Ar 2 ~ n Co _1 increases faster than the first one. Therefore the solution describing the 
space expansion can be possible only for A < 0. In this case at t — > oo, for m = as 
well as for m ^ 0, one can get the asymptotic representation (3.17) of the solution. 
This solution, as for the choice 1 < n < 2, provides asymptotically isotropic expansion 
of the B-I space. For t — > in this case we shall get only singular solution of the form 
(3.19). 

II. We study the system when Ljv = F( J), which means in the case considered D = 0. 
Let us note that, in the unified nonlinear spinor theory of Heisenberg the massive term 
remains absent, as according to Heisenberg, the particle mass should be obtained as a 



result of quantization of spinor prematter [24]. In nonlinear generalization of classical 



field equations, the massive term does not possess the significance that it possesses in 
linear one, as by no means it defines total energy (or mass) of nonlinear field system. 
Thus without losing the generality we can consider massless spinor field putting m = 
that leads to $ = 0. This assumption metamorphoses (2.16) to get 

P{t) = D q /t, D = const. (3.20) 

The system of equations (2.19) in this case reads 

U[-U 3 = 0, Uz-U 4 = 0, 

U^ + U! = 0, U' 4 + U 2 = 0. (3.21) 

Differentiating the first equation of system (3.21) and taking into account the third 
one we get 

U'l + V x = 0, (3.22) 

which leads to the solution 

U x = D x e ia + iD 3 e~ ia , U 3 = iD x e i(J + D 3 e~ ia . (3.23) 
Analogically for U 2 and C/ 4 one gets 

U 2 = D 2 e ia + iD 4 e- ia , U 4 = iD 2 e ia + D 4 e~ ia , (3.24) 
where are the constants of integration. Finally, we can write 

V x = (l/v^XZ^ + ^se—), V 2 = (l/V^)(D 2 e™ + iD 4 e- i n, 

V 3 = {l/^){iD x e i(T + D 3 e~ ia ), V 4 = (l/Vr){iD 2 e ia + D A e~ i(T ). (3.25) 
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Putting (3.25) into the expressions (2.15) one finds 

P = 2 {D{ + D\ - D\ - DD/t. (3.26) 

Comparison of (3.20) with (3.26) gives D = 2 (Dj + D\ - D 2 - D\). 
Let us now estimate r using the equation 

f/r = 3/cA(n- l)P 2n , (3.27) 

where we chose L N = \P 2n . Putting the value of P into (3.20) and integrating one 
gets 

f 2 = -3K\D 2 n T 2 ~ 2n + y 2 , (3.28) 

where y 2 is the integration constant and can be defined from (2.27): y 2 = X/3 > 0. 
The solution to the equation (3.28) in quadrature reads 

<lT = t. (3.29) 



y/ ' -3K\Dl n T 2 - 2n + y 2 

Let us now analyze the solution obtained here. As one can see the case n — 1 is the 
linear one. In case of A < for n > 1 i.e. 2 — 2n < 0, we get 



r(t) | t ^o~ [( V^IAI/^] 1 /™, 

and 

r | t ^oo~ \/ 3/q/ 2 1, 

it means that for the term L^v considered with A < and n > 1 the solution is initially 
singular and the space-time is anisotropic at t — > oo. Let us now study it for n < 1. In 
this case we obtain 



l/n 



t It— >o~ \/3Ky 2 t, 

and 



r t- 



[(V^jA|I>?n)t] 



The solution is initially singular as in previous case, but as far as l/n > 1, it provides 
asymptotically isotropic expansion of B-I space-time. 

III. In this case we study = F(I, J). Choosing 

L N = F(K±), K+ = I + J = I V = -I A , K. = I-J = I T , (3.30) 
in case of massless NLSF we find 

V = 2SF K± , g = ±2PF K±) F K± = dF/dK±. 
Putting them into (2.16) we find 

S' 2 ±P 2 = J D ± . (3.31) 
Choosing F = \K± from (2.21) we get 

f = 3k\(ti - 1) Dl T 1 " 2 ", (3.32) 
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with the solution 

T 



n ~ x dT 



(3 33) 

I g 2 T 2n ~ 2 - ZkXDI 

where g 2 = X/3. Let us study the case with A < 0. For n < 1 from (3.33) one gets 

r{t) |^o~ Qt -> 0, (3.34) 
i.e. the solutions are initially singular, and 

r{t) [^IMDl)^, (3.35) 

which means that the anisotropy disappears as the Universe expands. In case of n > 1 
we get 

r(t) | t _*o« -> 0, 

and 

i.e. the solutions are initially singular and the metric functions a(t) , b(t) , c(t) are dif- 
ferent at t — > oo, i.e. the isotropization process remains absent. For A > we get 
the solutions those are initially regular, but it violates the dominant energy condition 
in Hawking-Penrose theorem f[8| . Note that one comes to the analogical conclusion 
choosing L N = XS 2n P 2n . 



4 Analysis of the results obtained when the B-I 
Universe is filled with perfect fluid 

Let us now analyze the system filled with perfect fluid. Let us recall that the energy- 
momentum tensor of perfect fluid is 

T; {m) = (p + e)uX - by = (e, -p, -p, -p). (4.1) 

As we saw earlier the introduction of perfect fluid does not change the field equations, 
thus leaving the solutions to the NLSF equations externally unchanged. Changes in 
the solutions performed by perfect fluid carried out through Einstein equations, namely 
through r. So, let us first see how the quantities e and p connected with r. In doing 
this we use the well-known equality T^. u = 0, that leads to 

j t (re) + fp = 0, (4.2) 

with the solution 

lnr = - / (4.3) 



Recalling the equation of state p = £e, < £ < 1 finally we get 



T° - f - T 1 - T 2 - T 3 - -7i - °^ (A A) 

- L 0(m) — fc — i+£> l(m) — - t 2(m) — -^(m) — — T l+? ' l^-^y 1 
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where e is the integration constant. Putting them into (2.21) we get 



t _ 3k (£ - l)e 



T 



(s+ i, < 4 ' 5 > 



which shows that for stiff matter (£ = 1) the contribution of fluid to the solution is 
missing. Let us now study the system with nonlinearity type I. In this case we get 

/ dT = ±y/3ih. (4.6) 

J JmC G T-\C^/T^- 2 ) +e Q r^) +g 2 



As one can see in case of dust (£ = 0) the fluid term can be combined with the massive 
one, whereas in case of stiff matter (£ = 1) it mixes up with the constant. Analyzing 
the equation (4.6) one comes to the conclusion that the presence of perfect fluid does 
not influence the result obtained earlier for the nonlinear term type I. One comes to the 
same conclusion analyzing the system with perfect fluid for the other types of nonlinear 
terms considered here. At least both at t — » and at t — > oo the key role is played by 
the other terms rather than the term presenting fluid. 



5 Conclusions 

Exact solutions to the NLSF equations have been obtained for the nonlinear terms being 
arbitrary functions of the invariant I = S 2 and J = P 2 , where S = tptjj and P = iip^y 5 ^ 
are the real bilinear forms of spinor field, for B-I space-time. Equations with power 
nonlinearity in spinor field Lagrangian = XS n , where A is the coupling constant, 
have been thoroughly studied. In this case it is shown that equations mentioned possess 
solutions both regular and singular at the initial moment of time for n > 2 . Singularity 
remains absent for the case of field system with broken dominant energy condition. It 
is also shown that if in the NLSF equation the massive parameter m ^ and n > 2 
then at t — > oo isotropization of B-I space-time expansion takes place, while for m = 
the expansion is anisotropic. Properties of solutions to the spinor field equation for 
1 < n < 2 and < n < 1 we also studied. It was found that in these cases there 
does not exist solution that is initially regular. At t — > oo the isotropization process of 
B-I space-time takes place both for m ^ and for m — 0. In case of nonlinear term 
Ln = XP 2n , we found the solutions those are initially singular and the isotropization 
process of B-I space-time depends on the choice of n. For = A(/± J) n we obtained 
the solutions those may be initially singular or regular depends on the sign of coupling 
constant A, but isotropization process depends on the value of power n. It is also shown 
that the results remain unchanged even in the case when the B-I space-time is filled 
with perfect fluid. 
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